PbTe/PbSnTe heterostructures as analogues of topological insulators 
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We investigate theoretically the PbTe/Pbi-^Sn^Te heterostructure grown in [111] direction, 
specifically a quantum wall (potential step of width d) of PbTe embedded in Pbi-^Sn^Te. For 
x large enough to lead to band inversion, and for large d, there are well-known gapless interface 
states associated with four L valleys. We show that for d m 10 nm the three pairs of states from 
oblique valleys strongly couple, and become gapped with a gap ~ 10 meV. On the other hand, 
the interface states from the [111] valley are essentially uncoupled, and they retain their helical 
character, remaining analogous to states at surfaces of thin layers of three-dimensional topological 
insulators. 



Introduction. A large attention has been recently de- 
voted to topological insulators (TIs), materials in which 
the topology of the bulk bandstructure guarantees the 
existence of robust (against non-magnetic disorder) spin- 
nondegenerate (helical) edge/surface states [HE]- TIs are 
essentially a subclass of narrow gap semiconductors, since 
their bandgaps must be smaller than typical energy scale 
of relativistic corrections to the bandstructure (<1 eV). 
In fact, the first two-dimensional TI was realized [3] in 
a quantum well of HgTe/HgCdTe, i.e. a heterostructure 
based on well-known compound narrow-gap semiconduc- 
tors. Helical edge modes have also been recently reported 
in a heterostructure of InAs/GaSb [3]. 

Part of interest in TIs stems from the fact that they 
represent a new topological phase of noninteracting elec- 
trons: the TI character of a material is its bulk property, 
nontrivially encoded in the wavefunctions of the occu- 
pied (valence band) states. However, it is the presence 
of the helical edge/surface states which leads to observ- 
able consequences. In three dimensional strong (weak) 
TIs [5j they consist of an odd (even) number of spin- 
nondegenerate Dirac cones, and many interesting effects 
were predicted to occur when these states become gapped 
due to a perturbation [TJ[2] (e.g. due to proximity to a su- 
perconductor [5] or a ferromagnet [7JI1])- By now many 
materials have been shown to be 3D TIs (see [T] and ref- 
erences therein), with Bi 2 Se3 and Bi 2 Te3 gathering the 
most attention. Despite a recent progress in gating of 
thin layers of Bi 2 Se 3 [HUH], it would be desirable to in- 
vestigate the helical surface states in a well-known semi- 
conductor system, for which the growth and nanostruc- 
ture processing are already mastered. One such platform 
is strained bulk HgTe [H]. In this article we propose 
to focus on heterostuctures of compound semiconductors 
from the IV- VI lead chalcogenide family, (Pb,Sn)Te [12]. 

The superlattices of PbTc/Pbi^Sn^Te (with x < 0.18) 
were thoroughly characterized [13] ■ Due to the huge 
dielectric constant of PbTe, mesoscopic structures of 
PbTe/PbEuTe exhibit very robust ballistic transport 
properties [TU [15]. Both Mn and Eu can be incorpo- 
rated into PbTe, and exchange coupling between their 



spins and the band carriers was studied [16] . Further- 
more, there exist heterostructures of PbTe and europium 
chalcogenides (EuSe,EuTe) [T7J [TH], which are magnetic 
semiconductors with rich phase diagrams. Finally, a su- 
perconducting proximity effect has been recently shown 
to exist at In-PbTe junctions [15] . All these features sug- 
gest that IV- VI based heterostructure, if it hosted gapless 
helical states, could be a good candidate for investiga- 
tions of phenomena related to TIs. 

When changing the Sn content in Pbi-^Sn^Te alloy, 
band inversion between the topmost valence band and 
the lowest conduction band occurs at x ~ 0.37. In an un- 
strained sample this happens simultaneously at the four 
L points at the edges of the Brillouin zone. Since the 
topological class of the bandstructure changes at the in- 
stances of local bandgap closing [20] , after an even num- 
ber of local band inversions this class remains the same. 
In fact, both PbTe and PbSnTe are in a topologically 
trivial class [5], although very recently it was argued 
[2"T] that SnTe is a topological crystalline insulator |22j . 
i.e. that it supports an even number of surface states for 
specific surface orientations. However, we focus here on 
the mutual band inversion between PbTe and PbSnTe, 
which leads to existence of interfacial helical states, pre- 
dicted theoretically 25 years ago [2"3Tf2"6"] . 

Such an interface, with four Dirac cones, is analogous 
to a surface of a weak TI. Recent research on weak TIs 
suggests that these cones are in fact protected against 
time-reversal perturbations as long as these are not peri- 
odic and commensurate with the lattice spacing [371 US] > 
a situation which is rather improbable at an epitaxial in- 
terface of two materials with the same lattice structure. 
It would however be even more interesting if one could 
obtain a single gapless helical state at an interface. As we 
show in this article, this can be brought by a confinement 
effect in a properly designed PbTe/Pbi^Sn^Te het- 
erostructure grown in [111] direction. This happens due 
to an anisotropic energy structure near the four L points, 
which leads to a much stronger coupling of oblique valley 
states from opposite interfaces, compared to the [111] val- 
ley along the growth direction (a similar coupling opens 
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up a gap in a thin layer of a TI [29-32 ]). With a proper 
choice of strain, composition, and layer width, we ob- 
tain practically gapless helical states for [111] valley, with 
Dirac points located in the finite-width-induced gap in 
the remaining three valleys. 

The Hamiltonian. The electronic states near the L- 
point extrema in Pbi-^Sn^Te are described by Dim- 
mock's Hamiltonian |12U13| . Since we focus on HS grown 
in [111] direction, we choose the z axis parallel to [111]. 
The Hamiltonian in the basis of 4| (X - iY) \\), i (X + 

*Y) \l),^(S x -iS Y ) It), ^(Sx+iSy) II) (where S x .y 
transform like X, Y but do not change sign under reflec- 
tion) is then given by 
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where e(k) = V (x) + D 1 k 2 + D 2 k\ , M (k) = A (x) + B 1 k 2 + 
Bik\ and Q = Po • k. For [111] valley Po is a ma- 
trix with (v±,vx,v\\) on the diagonal. We stress that a 
realistic large ratio of v±/vn is crucial for our considera- 
tions. For remaining oblique valleys we replace Po with 
P = SPoS -1 , where S is a transformation matrix between 
the coordinate systems associated with the oblique valley, 
{e a /} with e 2 < along the valley direction, and {e a } with 
e z along the growth axis: (e 

) = ( 

We take A(x)= 0.095 — 0.26a; eV, corresponding to the 
gap inversion at x c w 0.37, and V{x) = 0.125x eV [13] , 
We neglect the composition dependence of the remaining 
paramters which is, if any, very weak in the considered 
regime of x < 0.5 Q21 [13], and we take D l( ^ = 1.65 (9.25) 
eVA 2 , B 1{2) = 6.15(49.35) eVA 2 , U|| = 1.44 eVA, and 
v_l=4.76 eVA H3J. 

The strain described by tensor e is accounted for by 
H s = diagj^Tre + D c u e zz , D^Tie + D v u e zz }, in which 
the two entries correspond to 2 x 2 blocks. Note that 
Hs gives strain-dependent corrections to A and V. The 
acoustic deformation potentials are taken as composition- 
independent: D% = -\m eV, D v d = -2.23 eV, £>£ = 2.07 
eV, and D£ = 2.62 eV [T3J. The lattice constants of PbTe 
(SnTe) are taken as a = 6.454 A (6.313 A). 

The above Hamiltonian is written in a form which al- 
lows one to immediately see that after a basis reordering 
it becomes the same as the recently studied model kp 
Hamiltonian describing the states near the T point in 
the Bi2Se3 family of materials [33]. There, the TI phase 
occurs when A • B\ ^ < |33J, i.e. when the conduc- 
tion and valence bands at T point become inverted. As 
mentioned before, in unstrained Pbi-^Sn^Te the inver- 
sion occurs simultaneously at four L points, and the bulk 
phase remain topologically trivial. 

Helical nature of interface states. We first analyze 
interface-bound states (IBS) localized at two decoupled 
heterointerfaces of a very wide PbTc/Pbi-^Sn^Te het- 
erostructure grown in [111] direction. We assume that 



for \z\ < d/2 we have the values of A > and V corre- 
sponding to PbTe, while A' < and V correspond to 
Pbo.54Sno.46Te. From the point of view of band lineup, 
this structure is in fact a finite width potential step - a 
quantum wall (QWa) instead of a quantum well (QWe) . 

We use the effective mass approximation and replace 
k z by —id/dz in Eq. ([I]). As long as the intervalley cou- 
pling (due to alloy disorder or atomic-scale interface re- 
construction) can be neglected, the presence of mutual 
band inversion at the interfaces guarantees the existence 
of helical Dirac- like IBS [531 [2H] • The dispersion and the 
spin structure of IBS is insensitive to exact profile of the 
interface (23J [2B], with only the z dependence of wave- 
functions being affected in the region of interdiffusion. 
For simplicity we assume an abrupt interface [24l EH E] , 
and use the boundary conditions of continuity of the 
wavefunction and its derivative (the latter only when we 
consider nonzero k 2 terms). 

Since our numerical calculations show that for real- 
istic parameters the k 2 terms are of minor quantita- 
tive importance, we present much simpler analytical re- 
sults obtained by neglecting them. We also focus on re- 
sults for [111] valley, which we label with superscript a 
(with f3 denoting the oblique- valley results). The en- 
ergies of the IBS are E±(k±) — Eq ± v±k± A V A , , where 



V 2 = (A - A') 2 - (V - V') 2 and E 
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valleys the dispersion is of course anisotropic in kj_). The 
eigenvectors are 4-spinors ip± (k_i_), with L/ R denoting 
the left (z — —d/2) and the right (z = d/2) interfaces, 
multiplied by appropriate functions / ± (z) describing 
the exponential decay away from a given interface. The 
inverse decay lengths are given by 
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in the middle layer and x'zp —X±(^ A') outside, with 
rj= (V — V')/(A — A'). Let us note that x± acquire an 
imaginary part at finite kj^ - however, what is crucial is 
the fact that at k± = 0, Xo 1S given by Eq. with v\\ 

replaced by | • / 8zr[ + v 2 



The two-dimensional subspaccs of IBS at the two inter- 
faces are spanned by 4-spinors {(f>^ R , <j)^ R }, for which 
the subscripts signify that their components correspond 
to well-defined spin states, although they remain lin- 
ear combinations of different orbitals. Defining 7 = 
arctan (1 + 77) /(l — if) we can write them as: ^ L ' R — 



[=Fz cos 7,0, sin 7, 0] J , (j)^ R = [0, ±icoS7, 0, smj] 1 . Us- 
ing these basis states we can write the 4-spinor parts of 
the IBS solutions: $£(kj_) = (±i<j>^ + e ie $f)/V2, and 

^(■k ± ) = (j>n±ie w 4 ) R )/V2, where e ie = {k x + ik y )/k ± . 

The helical nature of these states can be seen 
by calculating the expectation values of spin opera- 
tors: (V>± a x = =F cos 27 sin 6 , a v t/S^\ = 
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FIG. 1: (Color online) Solid line is the dispersion of the [111] 
valley interface-bound states for PbTe/Pbo.54Sno.46Te quan- 
tum wall grown in [111] direction, with PbTe widths d of 12 
and 100 nm (undistinguishable), plotted as a function of in- 
plane k vector. The dashed lines are the bulk band edges for 
Pbo.54Sno.46Te (the band edges of PbTe are outside of the 
shown energy region). There is a gap of 0.2 meV opened at 
the Dirac point for d=12 nm, see Fig. [3] 

± cos27COS#, i.e. the spin vector is perpendicular to k^, 
pointing clockwise (anticlockwise) for negative (positive) 
energy branch. The pattern on the R interface is re- 
versed. More generally |35| we can write an effective 
surface Hamiltonian in the {(f>^ R , (j>± } basis: 

HsJr? = ±v± A _ A , (a x k y ~a v k x ) , (3) 

with the a matrices operating in the respective two- 
dimensional spaces. Let us remark that the spin pattern 
for states from the oblique valleys is more complicated, 
with nonzero out-of-plane (a z ) polarization. 

All these results show complete analogy with the sur- 
face states in TIs such as Bi2Se3 [33]. The difference is 
the existence of four pairs of IBS, with the states from 
all the valleys coexisting at a given energy in the bulk 
gap, showing that the PbTe/Pbi_ x Sn x Te heterointerface 
is analogous to a surface of a weak TI. 

The effect of finite QWa width. We assume that the 
lattice mismatch between PbTe and Pbo.54Sno.46Te is 
shared in 4 : 1 ratio between the two materials. This re- 
quires an engineering of strain distribution by, e.g., choice 
of a substrate. If we neglect the k 2 terms the calculation 
is fairly straightforward [Ml 123 GH] ■ With the k 2 terms 
and for a general valley direction we have to resort to 
solving the problem numerically. 

In Fig.[l]we show the energy dispersion of the [111] val- 
ley states formed in the QWa. For both widths of d=12 
and 100 nm the IBS states from the two interfaces remain 
practically uncoupled: the gap opened for d = 12 nm is 
only «0.2 meV. On the other hand, as we show in Fig. [2] 
with decreasing d a sizable gap opens for the IBS-derived 
states in the oblique valleys, and it is ~ 10 meV for d=12 
nm. As shown in Fig. [3k, the gaps Eg ^ are visibly differ- 
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FIG. 2: (Color online) Solid line is the dispersion of the [111] 
valley interface-bound states for PbTe/Pbo.54Sno.46Te quan- 
tum wall with PbTe width d of 100 nm (upper panel) and 
12 nm (lower panel). The dashed (dotted) lines are the bulk 
band edges for Pbo.54Sno.4eTe (PbTe). The gap of ~ 10 meV 
opens for d=12 nm. 



ent. Thus, for most of practical purposes (e.g. for exper- 
iments done at temperature ~ 1 K) , in a structure with 
d « 10 nm the dispersion of the states from the [111] 
valley can be considered gapless, while the states from 
the oblique valleys are gapped. Let us also stress that at 
energies |_E±(kjJ — E Q \ ^> E gi the eigenstates are essen- 
tially the L/R interface IBS (the description in terms of 
massive Dirac fermions [31] is relevant only at small kj^ ) . 
Furthermore, as shown in Fig.[3]3, for the chosen QWa the 
[111] valley states are in the gap of the oblique valley so- 
lutions, so that disorder with typical energy scale smaller 
than this gap cannot mix the states from [111] valley with 
the others. Thus, for epitaxial interface of good quality 
we expect that the above-described system is basically 
equivalent to a layer of material such as Bi2Te3, with the 
only difference being the smaller size of the gap within 
which the surface states exist (~ 10 meV vs ~ 0.3 eV). 

The above results can be understood using a calcula- 
tion with the k 2 terms neglected and d assumed to be 
large. We orthogonalize the IBS and form their symmet- 
ric and antisymmetric combinations, approximating the 
exact solutions in large d limit. We calculate the gap by 
evaluating the matrix elements of the QWa Hamiltonian 
between them. The d-dependence of Eg ^ comes from ex- 
ponential decay of the IBS in the step region, described 
by inverse lengths \o (calculated for k± = 0) given in 



4 



100 




10 20 30 40 50 

d [nm] 







(b) 




\oblique valleys 






\mi]valie>\___^ 










) 


20 


40 




d [nm] 





FIG. 3: (Color online) (a) The energy gaps of the IBS- 
derived states from [111] valley and the three oblique valleys. 
The solid lines are the results of an exact calculation, while 
the dashed lines are obtained with Eq. Q, which neglects 
the k 2 terms in the Hamiltonian. (b) Band edge energies for 
IBS-derived states. The Dirac point of [111] valley states is 
located in the gap of the states from the remaining valleys. 

Eq. ([5]) and below. In this way we obtain the bandgaps 
(see also [24"] ) 

E */e = 43 ? _g^_ eM-x a/0 d). (4) 

The predictions of this formula agree well with the results 
of the full calculation shown in Fig. [3^,. 

Discussion. The above choice of a heterostructure was 
dictated by strain dependence of energies of band ex- 
trema in Pbi^^Sn^Te. As shown in Fig. |4j when the 
lattice constant of Pbi-^Sn^Te is adjusted to the PbTe 
constant in a heterostructure grown in [111] direction, 
with the realistic values of deformation potentials [13] we 
obtain that the gaps at [111] L point and the remaining 
ones. One can see then that in the Pbi-^Sn^Te quan- 
tum well (with PbTe being the barrier), it is practically 
impossible (without applying a compressive strain) to ob- 
tain a Dirac point of [111] valley IBS within the gap 
of the oblique valley IBS-derived states. For the QWa 
with material parameters adopted here, we only require 
a modest amount of additional tensile strain, which could 
be provided at low temperatures by BaF 2 substrate typ- 
ically used for IV- VI compounds [13]. Additional tun- 
ing of band edge offsets can also be obtained by using 
quanternary alloys and replacing a fraction of Te by Se. 
Note however that if band-inverted Pbi-^Sn^Te itself 
possesses gapless surface states [2TJ, then the choice of 
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FIG. 4: (Color online) The energies of extrema of conduc- 
tion and valence bands for Pbi-^Sn^Te lattice matched to 
PbTe (grown in [111] direction), calculated using Eq. |TJ as a 
function of x. 



a compressively strained Pbi-^Sn^Te /PbTe quantum 
well will be more experimentally practical. 

Let us note that a result qualitatively the same as 
shown in Fig. [4] can be obtained for Pbi-^Sn^Te uni- 
axially strained in [111] direction. The non-overlapping 
of bandgaps in the range of x corresponding the inverted 
oblique valley and non- inverted [111] valley implies that 
applying such a strain, as proposed in [5] , cannot change 
Pbi-^Sn^Te into a strong TI, but into a topological 
semimetal. 

Another possible practical obstacle in investigation of 
helical states in the proposed structure is the issue of 
p-type self-doping of Pbi-^Sn^Te with x > 0.2. This 
is analogous to the situation in materials from Bi2Se3 
family, which possess large bulk conductivity when as- 
grown, see e.g. [35]. At x = 0.25 holes were compensated 
by indium doping [37], but it remain to be investigated 
how effective this would be at x~0.4 considered here. 

Conclusions. We have calculated the bandstructure 
of a heterostructure based on PbTe and band-inverted 
Pbi^Sn-rTe, in which PbTe forms a quantum wall 
(QWa) of finite thickness d. The "normal" quantum 
well states are then absent, and the localized solutions 
appear due to the mutual band inversion between the 
constituent materials. For very large d these solutions 
correspond to four pairs of Dirac cones located at the 
two interfaces [23H25] . As the QWa grown in the [111] 
direction is narrowed, these states become gapped due 
to their overlap. We have found that the opened gap is 
at least two orders of magnitude larger for states asso- 
ciated with three oblique valleys compared to the states 
from the [111] valley. Thus, after properly choosing the 
QWa composition, its width, and the strain distribution 
in it, one can obtain a practically gapless pair of helical 
Dirac states from the [111] valley, with their Dirac point 
located in the gap of the states from the other valleys. 
The IV- VI heterostructure can then serve as an analogue 
of a thin layer of strong topological insulator. 
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